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ABSTRACT 

We consider the application of high-pass Fourier filters to remove periodic systematic fluctuations from full-sky survey 
CMB datasets. We compare the filter performance with destriping codes commonly used to remove the effect of residual 
1/f noise from timelines. As a realistic working case, we use simulations of the typical Planck scanning strategy 
and Planck Low Frequency Instrument noise performance, with spurious periodic fluctuations that mimic a typical 
thermal disturbance. We show that the application of Fourier high-pass filters in chunks always requires subsequent 
normalisation of induced offsets by means of destriping. For a complex signal containing all the astrophysical and 
instrumental components, the result obtained by applying filter and destriping in series is comparable to the result 
obtained by destriping only, which makes the usefulness of Fourier filters questionable for removing this kind of effects. 
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1. Introduction 

Controlling systematic errors is crucial to CMB anisotropy 
measurements, and will likely become one of the most im- 
portant factors limiting the accuracy of future polarisation 
experi ments performed with ultra-high sensitivity detector 



arrays ( Mennella et al.|["2004| ). Second and third-generation 
CMB space missions (namely WMAFQ and P/anc/Q have 
been designed with tight systematic error control require- 
ments calling for the development of spacecraft with stable 
thermal interfaces, optimized orbits, and scanning strate- 
gies and, in the case of Planck, advanced cryogenic sys- 
tems for cooling high sensitivity differential receivers and 
bolometers with a high degree of stability. The sensitiv- 
ity achievable by present and future experiments on CMB 
anisotropy and polarisation requires that residual system- 
atics be controlled at the sub-/iK level. 

Cryo-coolers, in particular, can be the source of peri- 
odic systematic effects caused by fluctuations in the physi- 
cal temperature of the receivers and the warm electronics. 
Even after optimal in-hardware stabilization, a small level 
of residual fluctuations are generally present in the scientific 
data, with an impact that needs to be assessed to decide 
whether in-software removal should be applied before sci- 
ence exploitation. 

Because the thermal mass of both the satellite and the 
instrument damps high frequency temperature fluctuations, 
these effects are characterised by a frequency spectrum 
dominated by low frequencies, i.e. <C 1 Hz, that propa- 
gate to the measured CMB maps leaving a signature af- 



Send offprint requests to: F. de Gasperin 

1 http://map.gsfc.nasa.gov/ 

2 http://planck.esa.int 



ter being reduced by redundant measurements performed 
in each pixel. From the point of view of scientific data 
quality, the amplitude of the angular power spectrum of 
this residual must be significantly smaller than the largest 
residual caused by 1/f noise fluctuations remaining after 
map-making. 1/f residual is present even when using opti - 
mal map-making approaches (see, e.g. 
and appear as correlated structure on 



level of white noise ( |Maino et al.||2002 ). 



Ashdown et al.|2007 ) 

arge scales below the 



Time-ordered data (TOD) from CMB experiments are 
generally processed before map-making to remove or reduce 
the contaminati on from spurious effects. Some codes, such 



as destriping ([Kurki-Suonio et al. 2009| Keihanen et al. 



200 4] IPoutanen et al.|2004[) an d TourleTtilters ( |Tristraml£ 
Ganga 2007; |Hivon et al.|2QQ2| ) are able to achieve this with 
no or minimal assumptions about the effect to be removed. 
They are robust, relatively easy to implement, and widely 
used as standard tools in CMB data analysis. Other codes 
make strong assumptions about the effect and use other 
data (e.g. from temperature and/or electrical sensors) to 
detect and remove the spurious signals. These can provide 
excellent results provided that accurate complementary in- 
formation (e.g. from house-keeping telemetry such us the 
temperature recorded by a sensor on the instrument focal 
plane) is available. To remove thermal effects, for exam- 
ple, an effective use of temperature sensor data in non- 
blind codes calls for detailed knowledge about the ampli- 
tude and phase of thermal damping between the position 
of the temperature sensor s and the detectors (see, for ex- 
ample, Tomasi et al.|20fO ). 

In this paper, we consider the application of Fourier 
filters to CMB datasets, with particular reference to full- 
sky surveys performed from space. Fourier filters are well- 
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suited to remove spurious effects with known spectral shape 
and have been applied to data acquired by sub-orbital mis- 
sions. In the Boomerang experiment, for example, scan- 
synchronous effects were removed by filtering data with a 
high-pass filter with a cut frequency about seven tim es the 
scan frequency flHivon et alT]|2QQ2| [Masi et al.||2QQ6[ ). The 
application of Fourier filters to large datasets such as those 
produced by space surveys, however, is not straightforward. 
Aggressive filtering would cut the sky signal on large angu- 
lar scales that represent an essential part of full-sky CMB 
surveys. Furthermore, the large data stream size requires 
filters to be applied in chunks, introducing offsets that sub- 
sequently require normalisation. 

The objective of this work is to analyse the applica- 
tion of high-pass filters to the removal of slow periodic ef- 
fects from full-sky surveys datasets and compare the per- 
formance with destriping codes commonly used to remove 
the effect of residual 1/f noise from timelines. Data used 
in this work consist of simulations representing the typi- 
cal Planc k scanning strategy and P/ancfc-LFPj nois e per- 
formance flMennella et al.|2010| |Meinhold et al.|2QQ9| ) with 
a periodic fluctuations representing a thermal perturbation 
in the 20 K stage. The spectrum of the fluctuations is based 
on laboratory tests on the 20 K cooler, while the amplitude 
is adjustable according to the need of the simulation. In 
any case, the fluctuations are not representative of the ac- 
tual stability measured in-flight. We stress therefore that 
our analysis does not provide an assessment of the thermal 
systematic errors in Planck, but a tool to evaluate the ef- 
fectiveness of Fourier filters in removing periodic spurious 
fluctuations in CMB data. 



2. Residual of periodic systematic effects in CMB 
maps 

We introduce now some general considerations about the 
propagation of periodic fluctuations from measured time 
streams to CMB maps. Further details about this analytical 
description and a comparison with simulations performed 
i n the context of Planck can be found in Menne lla et al.l 
([2002). 



In what follows, we assume that each pixel of a sky ring 
is measured N times before the optical axis is repointed 
by an angle # rep . If we consider a periodic fluctuation of 
general shape in the detected signal (ST), we can expand it 
as the Fourier series ST — Y^j=-oo A? exp(z27r/jt), where fj 
represents the various fluctuation frequency components. 

To estimate how the amplitude of each component is 
reduced by the measurement redundancy provided by the 
scanning strategy, we divide the frequency spectrum of 
the systematic effect into two regions with respect to the 
scanning frequency, / sca n- (i) the low frequency region, 
with fj < / S can, and (ii) the high frequency region, with 

fj ^ /scan- 

In the low frequency region, each harmonic of amplitude 
Aj will be damped by the measurement redundancy by a 
factor proportional to sin(7r/j// scan ). In the high frequency 
region, instead, we differentiate between scan-synchronous 
(i.e. fj = kf scan ) and scan-asynchronous (i.e. fj ^ k / scan ) 
fluctuations. The first type will not be damped, as they 
are practically indistinguishable from the sky measurement; 



the second ones, instead, will be damped by the number of 
measurements for that pixel, that is a factor of the order of 
TV x # pix /# rep where 

$pix is the pixel angular dimension. 
If we also consider the additional reduction provided 
by the application of removal algorithms to the TOD, and 
denote with Fj the additional damping of each frequency 
component, then the final peak-to-peak effect of a generic 
signal fluctuation ST on the map can be estimated as 



(<5TP"P) map 
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3. Map-making and removal approaches 

Map-making is the process that combines the satellite 
pointing information and the instrument TODs int o a map 
flWright|1996l |Tegmark|1997a| |Stompor et al.|2002[ ) . Simple 
maps can be obtained by simply phase-binning data be- 
longing to the same sky ring and then by averaging binned 
data that are observed in the same sky pixel. This method 
produces a raw map that is usually affected by spurious 
signatures caused by 1// noise, long-term drifts and slow 
periodic effects. More sophisticated met hods can be divided 
into optimal, least squares map-making (I Wright et al.l 19961 



^atoli et al.l 



2001; 



Tegmark||1997a|bl |Borrill et aLl|200T 
Dore et al. 2001D, and approximate methods, such as de- 
striping (e. g.|Delabrou ille 1998; [Burigana et al.|l999|^ ' 
et al.|1999l|Revenu et al.|2000|[Keihanen et al.|2QQ4| 
Suonio et al. 120091). 
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Optimal methods require the inversion of large matri- 
ces or the use of iterative algorithms to avoid direct matrix 
inversions. Destriping algorithms, however, do not yield op- 
timal map s but are much less demanding on both memory 



and CPU (Poutanen et al. 2004 



Stompor & White 2004 



Efstathiou 2005[ |2007p and have been succes sfully applied 
to CM B data from the Planck space mission ( jZacchei et al.| 
|2011[ ). In this work, we produced maps using the Planck- 
LFI destriping code, which is briefly described in the next 
section. 



3.1. Destriping 

The destriping technique has been developed to reduce the 
effect of 1// noise fluctuations. The main requirement is 
that the knee frequency, /knee, must not be much higher 
than the scan frequency, / sca n; previo us simulation works 
perfo rmed in the context of Planck (Kurki-Suonio et al. 
2009) have shown that destriping is effective for knee fre- 



Low Frequency Instrument 



quencies up to 0.1 Hz, i.e. up to about six times the spin 
frequency. 

According to this approach, the contribution from 1/f 
noise or other slow fluctuations is approximated as a con- 
stant over a certain length of time. The code then estimates 
these constant baselines and removes them from the TOD. 
The baseline length (that was fixed to 60 seconds in our sim- 
ulations) can be chosen to optimise the cleaning, according 
to the time-scale of the systematic effect. 
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To identify the optimal baselines, ring crossings (i.e. 
common pixels observed from different scan circles) are 
searched into every baseline-long set of data. We indicate 
with Tij and n^i the observed signal and the white noise 
level for the pixel in the z-th row (circle) and j-th column 
(sampling along the circle) in matrices T and n, respec- 
tively. Baselines Ai are then recovered by minimising the 
quantity 



crossings 



( T i,l - Tj,m)Y 



n i,l 



n- 



(2) 



where the sum is over all the pairs of pixels present in the 
two different baseline-long set of data. Baselines are recov- 
ered by solving a set of linear equations and then can be 
easily subtracted from the data. 

Because destriping has also pr oven effective in remo ving 
slow periodic systematic effects ( |Zacchei et al. 2011| ), we 
also compare its performance with that of the high-pass 
filter. 



3.2. TOD processing before map-making 

Although destriping can reduce the impact of low frequency 
modes, greater suppression may be achieved using addi- 
tional algorithms that process TODs prior to map-making 
and can operate in "blind" or "non-blind" modes. Blind 
codes make no or little assumptions about the physical 
and/or statistical properties of the signals to be removed 
and use temporal information in the data. When combined 
with map-making, the spatial information (i.e. the corre- 
sponding pixel position for an observed data-point on the 
final map) can also contribute to the final suppression level 
of the systematic effect. 

Non-blind codes use the available house-keeping (H/K) 
data, providing auxiliary information such as currents and 
voltages to the various electronic devices, temperature sen- 
sors data etc. A simple example of non-blind analysis is 
the correlation between instrument output and H/K data. 
These methods are implemented by correlating the relevant 
H/K data with the radiometric output also using, where 
available, transfer functions obtained from dedicated tests. 
Another non-blind approach that has been studied within 
the Planck collaboration is o ne based on neural network 
algorithms ( |Maris et al.||2QQ4[ ). In this case, the neural net- 
work is "trained" using simulated sensor data together with 
radiometer output obtained via analytical transfer func- 
tions. As usual, after a "training" period, the network is fed 
with real instrument output and H/K data. Both methods 
clearly require H/K data from many sensors placed very 
close to the detectors on the focal plane. However, this is 
not always possible because of constraints on the number 
and position of the sensors in the instrument focal plane, 
which often limit the effectiveness and reliability of non- 
blind approaches. 



4. Data and procedures 

We describe now the procedure we used to simulate the 
time-streams and temperature anisotropy maps containing 
astrophysical and instrumental signatures. We also describe 
the high-pass filter used to remove the low frequency com- 
ponents of the instrumental periodic fluctuations. 



The simulations were executed considering the Planck 
scanning strategy and, in particular, instrumental parame- 
ters typical of the Planck-LFl 30 GHz receivers. A template 
for the periodic spurious fluctuation was provided by the 
expected temperature fluctuations in the 20 K LFI focal 
plane unit. 

4.1. Scanning strategy and sampling 

All the simulations were performed using the so-called "cy- 
cloidal" scanning strategy. According to this scheme, the 
satellite orbits around the L2 Lagrangian point of the Sun- 
Earth system and spins around its axis (at 1 rpm in our 
simulations), which is pointed towards the Sun with the 
solar panels keeping the pay load in shade. The telescope 
points at an angle of 85° with respect to the spin axis and 
sweeps the sky in near- great circles. The spin axis of the 
satellite also performs a slow circular path around the anti- 
solar axis, which is equivalent to a cycloidal path as a func- 
tion of ecliptic longitude. 

With a repointing of 2' every 48 minutes, the entire sky 
is observed by all detectors in the focal plane in about 7 
months. A small hol^Jnear the south pole was filled during 
the data-reduction procedure with values AT = K. 

All data-streams were produced with a sampling fre- 
quency of 32.5 Hz, which corresponds to that of the 30 GHz 
Planck-LFl radiometers and assuming that no gaps are 
present in the data streams. 

4.2. Time streams 
4.2.1. Astrophysical signal 

The sky emission template was provided by the Planck 
Sky Model (PSM), which consists of a set of tools that 
simulate the whole sky emission in the frequency range 
30 — 1000 GHz. In our study, we used a sky signal com- 
posed of the CMB and the diffuse foreground emission. 

The CMB was generated starting from the temperature 
powe r spectrum and fitting five years WMAP temperature 
data (Hinshaw et al. |2009[ ), while for the foreground emis- 
sions we assumed three galactic components: thermal dust, 
synchrotron, and free-free. The SZ effect and strong point 
sources were not included because their effect on the long- 
period variations are negligible. 

A main beam with #fwhm = 32.4 arcmin was convolved 
with the sky signal to obtain a realistic data-stream. The 
main beam was simulated with the software GRASP using 
measured feed-horn beam patterns and a detailed model of 
the Planck telescope, without considering the beam far side 
lobes. For more details about the LFI feed- horn patterns 
and estimated be ams in the sky, the rea der can refer to 
|Villa et al.|fl2009D and |Sandri et al.| ( |2010 ). 



4.2.2. Instrumental noise 

The instrumental noise was simulated as a combination of 
white noise plus al// a component with a power spectrum 
of the form 



p(f) 



knee 



(3) 



226 pixels in a map with N P u 



786432, that is 0.03%. 
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where a 2 is the white noise variance, /knee is the frequency 
at which white noise and 1/ f a contribute equally in power 
and a is the exponent in the power law that is usually in 
the range 0.5 < a < 2. 

The presence of l/f a noise in CMB data acquired 
by coherent receivers such as those used in Planck-LFI 
flBersanelli et al.||201Q[ ) and WMAP ( jBennett etHl|2003 ) 
is caused by gain and noise temperature fluctuations in the 
high electron mobility transistor (HEMT) amplifiers. These 
fluctuations can be highly suppressed in hardware b y adopt- 
ing efficient differential receiver schemes (see, e.g.: |Seiffert| 
|etaT1[2QQ2| [Mennella et al.||2QQ3| . 

Noise time-streams were simulated using the inverse 
Fourier transforming equation Q with the noise param- 
~ In particular, the noise spec- 



eters reported in Table 
tral density £Ti sec _ANT is related to the white noise vari- 
ance <j 2 by the equation a = 5Ti sec _ ant / V^nt , where 
Tint = V/samp * s ^ ne sample integration time. 

Table 1: Typical noise parameters for Planck-LFI 30 GHz 
receivers 



/knee 


a 


^Ti S ec-ANT 


/samp 


0.008 Hz 


0.69 


2.629 • 10" 4 Ks 1/2 


32.5 Hz 



4.2.3. Periodic systematic effect 

Periodic, "slow" spurious fluctuations in cryogenic mi- 
crowave receivers are often caused by temperature insta- 
bilities in the cryogenic system that propagate through the 
mechanical structure and couple with the measured signal. 
In coherent receivers, for example, this coupling is provided 
by the amplifier gain and noise temperature that oscillate 
synchronously with the environment physical temperature 
flTerenzi et aL|2009| ). 

In our work, we simulated a realistic periodic effect 
starting from th e expected stability behaviour of the Planc k 
sorption cooler ( |Bhandariet al. 2004 Morgante et al. 2009), 
which is the main cryogenic system on board the Planck 
satellite providing a 20 K stage for the LFI receivers and 
a 18 K precooling stage for the HFI 4 K cooler. The main 
frequencies of the temperature fluctuations depend on the 
(programmable) absorption-desorption period of both each 
single compressor and the whole 6-compressor assembly. We 
chose a typical value of 940 s for the single compressor pe- 
riod which leads to 5640 s for the main period of the whole 
assembly. The design temperature stability at the cooler 
cold end is ~ 0.1 K peak-to-peak. 

Details of the thermal damping provided by the instru- 
ment structure and of the radiometric coupling functions 
are beyond the scope of this paper a nd not discussed here . 
Interested readers can find details in Terenzi et al. (2009) 



and |Tomasi et al.| ( |2010[ ). 

In the right panel of Fig. [I] we show an example of the 
periodic effect in the receiver output in the time and fre- 
quency domains. The peak-to-peak effect in antenna tem- 



5 These parameters are representative of Planck-LFI 30 GHz 
receivers noise performa nces as measured during the instrument 



perature is ~ 4 mK and the two yellow lines in the spec- 
trum highlight the main frequency peaks of the systematic 
effect. In Fig. |2j the same effect is shown when projected 
onto a map. The visible stripes are the signature of the 
systematic effect that has been damped by the redundancy 
in the scanning strategy (see Sect. |2| to a peak-to-peak of 
< 70 //K. 

4.3. High- pass filter 

The power spectra of both the astrophysical signal and the 
periodic effect shown in Fig. [l] highlight their different dis- 
tribution in the frequency domain. Because of the scan- 
ning strategy, the astrophysical signal lies mainly at the 
frequency of 1/60 Hz and overtones; the periodic effect, in- 
stead, is characterised by fluctuations at lower frequencies. 
A high-pass filter appears a natural solution to remove the 
systematic effect without affecting the astrophysical com- 
ponent. 

The filter action H(f) in the frequency domain, can be 
described by the equation 

[£(/)]hpf = sky + [S(f)]p + [S(f)]l/f ), (4) 

where [£(/)] hpf is the spectrum of the filtered signal, 
[S(f)] s ky, [*§(/)] p > and [*§(/)] i// are the Fourier transforms 
of the astrophysical component, the periodic systematic ef- 
fect, and the instrument noise. 

The high-pass filter used in our w ork is a simple real 
Butterworth filter ( |Butterworth| p"930 ) of the form 

(///cut) n 



H(f) 



(// /cut)^ 



(5) 



groun d test campaign (Mennella et al. 2010 Meinhold et al. 
2009) 



where / cut represents the filter cut-off frequency and n the 
slope. This filter removes almost all the signal at frequencies 
/ < /cut, leaving the complementary range of frequencies 
nearly untouched. Because the Fourier transform is linear, 
the filter can be applied to the different signals separately 
to remove the systematic effect and leave the astrophysical 
signal unchanged. 

A C++ code was developed and run on a 32-node 
cluster of Intel Xeon 3.00 GHz processors with 2.00 GB 
SDRAM each; on this system, the code can process a seven- 
month dataset for a single detector in few minutes. 

Our first study has been a sensitivity analysis to deter- 
mine the optimal cut frequency, / cut , and the filter step, 
TV, i.e. the length of data that the filter will process at each 
step. This last parameter, in particular, is required to avoid 
cutting the sky signal on large angular scales. Moreover, 
several months of data cannot be filtered in a single step 
because of memory limitations. The filter slope, n, was fixed 
in all our simulations at n = 24, a value that makes the fil- 
ter very steep yet avoids ringing side lobes that would be 
caused by a step function. 

In the sensitivity analysis, we optimised each parameter 
independently, while keeping the other two fixed at a ref- 
erence value. We first generated a data-stream with all the 
components described in Sect. |4.2| and with a periodic effect 
amplitude 100 times larger than the one shown in Fig. [I] 
We then calculated residual maps by subtracting the astro- 
physical signal map from the cleaned all-components map, 
i.e. 



Mt 



(residual) 



= [M ( 



sky+systematics) J HPF 



M { 



(sky) j 



(6) 
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Time [s] 




Freq [Hz] Freq [Hz] 

Fig. 1: The astrophysical signal (left) and the periodic systematic effect (right) in the simulated data-streams. The top 
panels of the two figures present the signals in the time domain, while the bottom panels show their power spectra. The 
green line is located at the frequency of 1/60 Hz, while the two yellow lines indicate the two main frequencies of the 
periodic effect. The red dashed line shows the shape of the high-pass filter; the filter curve is not to scale, its true value 
on the y-axis ranges from to 1. 




-3.8e-05 i 



2.8e-05 [K] 



Fig. 2: Map of the periodic systematic effect for a 30 GHz radiometer after a survey of 7 months showing structures 
caused by the slow signal periodic fluctuations. 



where M( sky ) is the map containing only the sky signal, 
while [M (sky+systematics) ] HpF is the map of all the compo- 
nents (i.e. astrophysical, instrument noise and the periodic 
effect) after high-pass filtering. We finally compared the an- 
gular power spectra, C^ es , calculate^ from residual maps 
obtained with different values of the considered filter pa- 
rameter. 



6 The angular power spectra were co mputed using the 
HEALPix ( jhttp : //healpix . jpl . nasa . gov)] utility anafast. 



In Fig.[3j we plot C\ es for various values of the multipole 
£ as a function of the cut-off frequency / cut . The n and N 
parameters were fixed, in these runs, at the values of n = 24 
and N = 24 hours, respectively. Smaller residuals clearly 
imply that the filter has a smaller effect on the astrophysical 
signal and that there is a larger systematic effect removal. 

If the cut-off frequency is too low, then the filter is in- 
effective in removing the periodic effect components at fre- 
quencies ~ 0.17 mHz (= 1/5640 s _1 ) and ~ 1 mHz 
(= 1/940 s _1 ), as reflected in the increase of C\ es at low 
frequencies 
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Fig. 3: Residual C\ es at multipoles ranging from 1 (red) 
to 500 (blue) as a function of the cut-off frequency / cut . 
Vertical dotted lines represent the values of / cut that have 
been tested and the black dashed line is the optimal one. 
The horizontal dotted line represents the mean of the 500 
multipoles. In these runs, we have multiplied the amplitude 
of the periodic systematic effect by a factor 100 to enhance 
its effect. 



Fig. 4: Residual C\ es at various multipoles (ranging from 1 
to 600 and averaged in groups of 40) as a function of the 
filter step, N. Red curves represent the residual after fil- 
tering a data- stream containing only the systematic effect, 
while blue curves show the residual obtained after filter- 
ing the astrophysical signal. Vertical dotted lines represent 
the tested values and black dashed line is the selected one. 
Horizontal coloured dotted lines are the mean of all 600 
multipoles and the black line represents the sum of the red 
and blue dotted lines. 



When / cut approaches the scan frequency, the filter then 
causes a significant removal of the first sky-signal peak, 
located at 0.016 Hz (= 1/60 s" 1 ). This results in Cf s 
increasing at values of / cut > 0.016 Hz. For even higher 
values of / cu t? the filter starts to remove also overtones, 
causing a second increase in C\ es . An optimal value for the 
cut frequency has been fixed at / cut = 0.007 Hz. 

The effect of the sensitivity analysis on the filter step is 
shown in Fig. [4] (see caption for details) . The results show 
that N must be chosen to be high enough to make the filter 
capable of detecting the systematic effect but low enough to 
avoid excessive removal of long-period astrophysical signals. 
In particular, we see a rather flat minimum in the residual 
effects up to a filter step of N = 24 hours, which is the 
value chosen for our simulations. 

Another degree of freedom that has been studied is the 
"baseline removal" i.e. the removal of the zero frequency 
component from the filtered spectrum. In this case, the 
choice also comes from a trade-off between an effective 
cleaning of the periodic effect (that benefits also from the 
baseline removal) and the requirement that the filter must 
not alter the astrophysical signal (which is affected by the 
baseline removal). The overall effect of the two different 
choices is discussed in Sect. 03 



5. Results and discussion 

We compare the Fourier filtering (with and without baseline 
removal) with the destriping process and a combination of 
the two, in removing the periodic spurious signal shown in 
the right panel of Fig. [I] 

We first compare results obtained with the two different 
approaches in two simple cases: (i) a data stream containing 
only the periodic effect and (ii) a data stream containing 



only the astrophysical signal. We then perform the same 
comparison by also taking into account the 1/f noise. 

5.1. Filter applied to periodic signal alone 

In Fig. [5j the periodic systematic effect is projected onto a 
HEALPix with N s ^ e = 256 (corresponding to a pixel size of 
~ 13.7 arcmin) before and after the application of various 
removal algorithms. The peak-to-peak effect on the map is 
w 66 /iK, that is reduced to ~ 11 /iK after filtering and 
destriping. 

If we now consider the results obtained by filtering with- 
out destriping (maps (c) and (e)) we can see that: 



when the baselines are not removed (Fig. 5c), the resid- 
ual map displays stripes caused by long-period oscilla- 
tions that are not removed by the filter. These stripes 
can be removed by a subsequent application of the de- 
striping code; 

when baselines are removed, the filtered-only map 



(Fig. 5e) is similar to the destriped-only map (Fig. 5b) 
Additional destriping after filtering (Fig. [5f|) does not 
yield any measurable improvement. 

To study the direct impact of the periodic signal on the 
CMB, we compare power spectra as done in the top panel 
of Fig. [6] From the figure we see that best performance 
is obtained with the combination of filter and destriping. 
With the filter-only, in particular, baseline removal plays 
a key role, as it leaves stripes that otherwise need to be 
subsequently removed by destriping. If we use the filter + 
destriping combination, however, baseline removal becomes 
irrelevant and the residual after the combined cleaning pro- 
cedure is of the order of ~ 70% of the residual obtained by 
applying destriping only. 
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(a) Periodic systematic effect 





(b) After destriping 




(c) After filtering 



(d) After destriping and filtering 





(e) After filtering (baselines removed) 





(f) After destriping and filtering (baselines removed) 



Fig. 5: Maps of the periodic systematic effect before and after filtering- alone (with and without baseline removal), 
destriping-alone, and filtering + destriping. Colours are rescaled in every map to highlight footprints of any system- 
atic effects. In the above table, we summarise the map r.m.s. values and the corresponding dumping factors. 



In the bottom panel of Fig. [6j we compare the CMB 
po wer spectrum obtain ed from the seven year WMAP best 
fit ( Jarosik et al.||20TT| ) including cosmic variance with the 
residual caused by the periodic effect. The figure shows that 
in all cases the effect is at least two decades below the 
CMB spectrum and that the combination of filtering and 
destriping is able to reduce the residual effect by another 
order of magnitude. Similar comparisons at small angular 
scales are provided in Fig. [8] 



5.2. Filter applied to astrophysical signal alone 

In the next step we applied the same algorithms to a data 
stream containing only a simulated astrophysical signal 
composed by CMB + galactic diffuse emission. 



In Fig. [7| we show the power spectra of residual maps 
obtained after applying of the various filtering approaches. 
Residual maps are defined as 



Mi 



(residual) 



[M ( 



sky)J H PF/Dest 



Mi 



(sky) 5 



(7) 



where M( sky ) is the reference sky map and [M( sky )j HPF y Dest 
is the same map after the application of the cleaning code 
(with any of the tested combinations). 

The curves in Fig. [7] tell us that any kind of algorithm 
applied to clean systematic effects is going to also alter 
the astrophysical signal. The high-pass filter, in particu- 
lar, leaves a residual especially when baselines are removed. 
This is probably due to the strong difference in the mean 
value of the various chunks of data because the signal is 
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FlUar (with baaellne removal) = ■ 
Fitter (without baseline removal) = 
Fitter+ Destriping (with baseline removal) = 
Flltsr+ Destriping (without baseline removal) = 
Periodic systematic effect = ■ 




„ 10 1 



u io u 



CMB (WMAP 7 year best fit) + cosmic variance 



Periodic systematic effect 
Periodic systematic effect after filtering 
Periodic systematic effect after filtering and destriping 




Multipoles (I) 

Fig. 6: Top: Power spectra the periodic systematic effect be- 
fore (red curve) and after removal. Green: destriping only. 
Filtering without baseline removal - Cyan: filtering only, 
yellow: filtering + destriping. Filtering with baseline re- 
moval - Blue: filtering only, orange: filtering + destriping. 
Note that the yellow, orange, and blue curves are super- 
imposed, meaning that the three cases produced the same 
results. Bottom: Impact of periodic fluctuations on large 
scales compared to the CMB and cosmic variance obtained 
from a best-fit solution of the seven year WMAP data. 




CMB+Galaxy = 
Destriping = 
Filter (with baseline removal) = 
Filter (without baseline removal) = 
Flrtar+Deatriplng (with baseline removal) = 
Fllter+ Destriping (without baseline removal) = 




400 

Multipoles ffi 



mostly concentrated along some rings (blue patches in the 
top panel of Fig. [7]) and, therefore, in just some data- 
chunks. When the mean values of these data-chunks are 
removed, strong marked stripes are generated in the resid- 
ual map. From these results, it is clear that any application 
of Fourier filters to real data requires additional baseline 
normalisation such as that performed by destriping codes. 



Fig. 7: Top: residual CMB+Galaxy map after filtering with 
baseline removal. Middle: residual CMB + Galaxy map after 
destriping. Bottom: power spectra of residual CMB+galaxy 
maps after destriping (green), filter-only without baseline 
removal (cyan), filter-only with baseline removal (blue), fil- 
ter without baseline removal + destriping (yellow), filter 
with baseline removal + destriping (orange, below yellow 
curve) . 



5.3. Filter applied to the combined signal 

The third step was to apply various cleaning procedures to 
data-streams containing the sky signal, the 1/f component, 
and the periodic systematic effect and then compare the re- 
sults. Since all the steps (filtering, destriping/mapmaking, 
and spectrum extraction) are linear, we can work separately 
on the different components without any loss of informa- 
tion. 

In Fig. [8j we plot the residual power spectra for ev- 
ery signal separated by cleaning method. Fig. [8a| shows the 
power spectra of the various signals when no cleaning is 



applied. In the other panels, we report the same spectra 
after applying various cleaning techniques together with 
the residual coming from the modification of the sky signal 
caused by the software removal. The green curve in Fig.|8b[ 
in particular, shows the residual level from 1/f noise after 
destriping. Our results indicate that both destriping alone 
and the combination of filter and destriping can suppress 
the effect of the periodic signal at a level that is at least 
two order of magnitudes below this residual. 

Fig. 8c (high-pass filter) again highlights the need to use 



a destriper in order to renormalise offsets generated by the 
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CMB + Galaxy = 
CMB+GaJaxy - Destrlplng = 
Periodic systematic effect - Destrlplng = 




(a) Signals 



(b) Destriper effect 




CMB + Galaxy = 

CMB+Qalaxy- Fllter+ Destrlplng {without baseline removal) = 
:lc effect - Hfter+ Destrlplng (without baseline removal) = 
1A-Filter+Destriping (without baseline removal) = 




(c) Filter effect 



(d) Filter + Destriper effect 




CMB + Galaxy = 

CMB+Qalaxy- Flrter+Deatrlplng (with baseline removal) = 
' ks effect - Filter + Destrlplng (with baseline removal) = 
1/f-Fifter+Destriplng (with baseline removal) = 




(e) Filter (no baselines) effect 



(f) Filter (no baselines) + Destriper effect 



Fig. 8: Power spectra main components before (a) and after cleaning with destriping (b), high-pass filter (c) and code- 
combination (d). The last line shows the effect when the baselines are removed by high-pass filter (e) and the code- 
combination (f). 



filter. A comparison of Figs. 8b (destriping) and 8d (code- 
combination) shows that the damping of the systematic 
effect obtained by applying destriping only is essentially 
the same as that obtained by combining the two codes. 



The two bottom plots show the effects of the high-pass 
filter and the code-combination with baselines removal. The 
filter-only baseline removal again produces a larger resid- 
ual, while no significant difference is seen when the filter is 
combined with destriping. 



6. Conclusions 



We have analysed the use of a high-pass Fourier filter to 
clean full-sky CMB datasets from periodic systematic ef- 
fects. As a test case, we have considered the baseline Planck 
satellite scanning strategy and instrument parameters typ- 
ical of the Planck-LFl 30 GHz receivers. A template of 
the spectral content of periodic spurious signal was derived 
from ground tests of the 20 K Planck-LFI focal plane unit. 
The effectiveness of Fourier filtering was then compared 
with destriping and with a third approach that combines 
filtering and destriping in sequence. 
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After a sensitivity study to define the optimal filter pa- 
rameters, we first processed a data stream containing only 
the systematic effect. The high-pass filter with baseline 
removal is effective in removing the periodic signal, with 
residual effects lower than those obtained by destriping. 

We then processed a data stream containing only the 
astrophysical signal; in this case the filter creates a new 
artefact, visible in the power spectrum at low multipoles 
and in the residual maps as sharp stripes. This residual ef- 
fect can be effectively removed by a subsequent application 
of destriping. 

The final p owe r spectrum (with all the components de- 
scribed in Sec. 4.2) obtained after filter + destriping differs 
from the one obtained after destriping by < 10 _16 K 2 for 
£ < 500 and < 10" 18 K 2 for £ > 500. This means that 
destriping is able to reproduce almost the same effect ob- 
tained by the filter + destriping combination. These results 
indicate that high-pass Fourier filters are not suitable for 
cleaning large CMB datasets from periodic systematic ef- 
fects, especially when destriping can be effectively applied. 



Poutanen, T., Maino, D., Kurki-Suonio, H., Keihanen, E., Sz Hivon, 

E. 2004, MNRAS, 353, 43 
Revenu, B., Kim, A., Ansari, R., et al. 2000, A&AS, 142, 499 
Sandri, M., Villa, F., Bersanelli, M., et al. 2010, A&A, 520, A7+ 
Seiffert, M., Mennella, A., Burigana, C, et al. 2002, A&A, 391, 1185 
Stompor, R., Balbi, A., Borrill, J. D., et al. 2002, Phys. Rev. D, 65, 

022003 

Stompor, R. & White, M. 2004, A&A, 419, 783 
Tegmark, M. 1997a, Phys. Rev. D, 56, 4514 
Tegmark, M. 1997b, ApJ, 480, L87+ 

Terenzi, L., Salmon, M., Colin, A., et al. 2009, JINST, 4, T12012 
Tomasi, M., Cappellini, B., Gregorio, A., et al. 2010, JINST, 5, T01002 
Tristram, M. & Ganga, K. 2007, Reports on Progress in Physics, 70, 
899 

Villa, F., D'Arcangelo, O., Pecora, M., et al. 2009, JINST, 4, T12004 
Wright, E. L. 1996, arXiv:astro-ph/9612006vl 

Wright, E. L., Hinshaw, G., & Bennett, C. L. 1996, ApJ, 458, L53+ 
Zacchei, A., Maino, D., Baccigalupi, C., et al. 2011, submitted to A&A 



Acknowledgements. The work in this paper has been supported by in 
the framework of the ASI-E2 phase of the Planck contract and has 
been carried out at the Planck-LFI Data Processing Centre located in 
Trieste at the Astronomical Observatory. Some of the results in this 
pape r have been derived using the HEALPix package ( |G6rski et al.| 



References 

Ashdown, M. A. J., Baccigalupi, C, Balbi, A., et al. 2007, A&A, 471, 
361 

Bennett, C. L., Bay, M., Halpern, M., et al. 2003, ApJ, 583, 1 
Bersanelli, M., Mandolesi, N., Butler, R. C, et al. 2010, A&A, 520, 
A4+ 

Bhandari, P., Prina, M., Bowman, R. C, et al. 2004, Cryogenics, 44, 
395 

Borrill, J., Ferreira, P. G., Jaffe, A. H., & Stompor, R. 2001, in Mining 
the Sky, ed. A. J. Banday, S. Zaroubi, & M. Bartelmann, 403 — h 

Burigana, C, Malaspina, M., Mandolesi, N., et al. 1999, arXiv:astro- 
ph/9906360 

Butterworth, S. 1930, Experimental wireless & the wireless engineer, 
536 

Delabrouille, J. 1998, A&AS, 127, 555 

Dore, O., Teyssier, R., Bouchet, F. R., Vibert, D., & Prunet, S. 2001, 

A&A, 374, 358 
Efstathiou, G. 2005, MNRAS, 356, 1549 
Efstathiou, G. 2007, MNRAS, 380, 1621 

Gorski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 
Hinshaw, G., Weiland, J. L., Hill, R. S., et al. 2009, ApJS, 180, 225 
Hivon, E., Gorski, K. M., Netterneld, C. B., et al. 2002, ApJ, 567, 2 
Jarosik, N., Bennett, C. L., Dunkley, J., et al. 2011, ApJS, 192, 14 
Keihanen, E., Kurki-Suonio, H., Poutanen, T., Maino, D., &; 

Burigana, C. 2004, A&A, 428, 287 
Kurki-Suonio, H., Keihanen, E., Keskitalo, R., et al. 2009, A&A, 506, 

1511 

Maino, D., Burigana, C, Gorski, K. M., Mandolesi, N., &; Bersanelli, 

M. 2002, A&A, 387, 356 
Maino, D., Burigana, C, Maltoni, M., et al. 1999, A&AS, 140, 383 
Maris, M., Mennella, A., Terenzi, L., et al. 2004, Memorie della Societa 

Astronomica Italiana Supplement, 5, 399 
Masi, S., Ade, P. A. R., Bock, J. J., et al. 2006, A&A, 458, 687 
Meinhold, P., Leonardi, R., Artal, E., et al. 2009, JINST, 4, T12009 
Mennella, A., Baccigalupi, C, Balbi, A., et al. 2004, Res. Devel. 

Astronomy & Astrophys., 2, 1 
Mennella, A., Bersanelli, M., Burigana, C, et al. 2002, A&A, 384, 736 
Mennella, A., Bersanelli, M., Butler, R. C, et al. 2010, A&A, 520, 

A5+ 

Mennella, A., Bersanelli, M., Seiffert, M., et al. 2003, A&A, 410, 1089 
Morgante, G., Pearson, D., Melot, F., et al. 2009, JINST, 4, T12016 
Natoli, P., de Gasperis, G., Gheller, C, & Vittorio, N. 2001, A&A, 
372, 346 



10 



